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ABSTRACT 

We present an analysis of the X-ray spectrum of the Local Bubble, obtained by simultaneously 
analyzing spectra from two XMM-Newton pointings on and off an absorbing filament in the Southern 
galactic hemisphere (b as —45°). We use the difference in the Galactic column density in these 
two directions to deduce the contributions of the unabsorbed foreground emission due to the Local 
Bubble, and the absorbed emission from the Galactic halo and the extragalactic background. We 
find the Local Bubble emission is consistent with emission from a plasma in collisional ionization 
equilibrium with a temperature log(TLe/K) = 6.06^ 04 and an emission measure J n^dl = 0.018 
cm~ 6 pc. Our measured temperature is in good agreement with values obtained from ROSAT All- 
Sky Survey data, but is lower than that measured by other recent XMM-Newton observations of the 
Local Bubble, which find log(7LB/K) as 6.2 (although for some of these observations it is possible 
that the foreground emission is contaminated by non-Local Bubble emission from Loop I) . The higher 
temperature observed towards other directions is inconsistent with our data, when combined with a 
FUSE measurement of the Galactic halo O vi intensity. This therefore suggests that the Local Bubble 
is thermally anisotropic. 

Our data are unable to rule out a non-equilibrium model in which the plasma is underionized. How- 
ever, an overionized recombining plasma model, while observationally acceptable for certain densities 
and temperatures, generally gives an implausibly young age for the Local Bubble (<6x 10 5 yr). 
Subject headings: Galaxy: general — Galaxy: halo — ISM: general — ISM: individual (Local Bubble) — 
X-rays: ISM 



1. INTRODUCTION 

The Local Bubble (LB) is a region of X-ray-emitting 
gas of ^100 pc extent in which the Solar System resides. 
Measurements of the interstellar Na I absorption towards 
456 nearby stars reveal that the Local Bubble resides in a 
cavity in the interstellar medium fISM; ISfeir et al.lll999f) . 
The idea of the Local Bubble originated in the late 1970s, 
in order to explain the observed anticorrelation between 
the intensity of the soft X-ray backgroun d (SXRB) and 
the Galactic hydrogen colu mn density Nn (jBowver et al.l 
[19681: iSanders et all 11971 . The data are inconsistent 
with the anticorrelation being due to absorption, as they 
require an interstellar absorpt ion cross-section th at is 
one-third of its expected value (|Bowver et al.lfl96ct ). and 
also dif ferent energy bands have the same d ependence 
on Nh (jSanders et alj|1977t Uuda et al.lll991l ). Instead, 
a so-called "displacement" model was proposed, in which 
the hot X-ray-emitting plasma (the Loca l Bubble) is in 
the foreground and displaces the cool gas ( Sander s et al.l 
Il977t iTanaka fc Bleekerlll977ft . In directions of higher X- 
ray intensity the Local Bubble is thought to be of greater 
extent, and so there is less cool gas (and hence lower Nh) 
in those directions. 

Numerous models have been proposed for the 
formation of the L ocal Bubb l e (for reviews, see 
iBreitschwerdl 119981: ICoxl 119981 : IBreitschwerdt fc Coxl 
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2004). These models may essentially be divided into 
two classes. In one class of model, the Local Bub- 
ble was carved out of the ambient ISM by a su- 
perno v a or series of supe r novae ( e.g. ICox fc Anderson! 
119821 llnnes fc Hartqulstl 119841: ISmith fc Coxl 120011 : 
lMafz-Apellanij|200H IBreitschwerdt fc de Avillezll2006h . 
The hot gas thus produced gives rise to the observed 
X-rays. If the last supernova was recent enough, one 
would expect the ions to be underionized. However, 
if the Local Bubble is old enough to have begun con- 
tracting, the ions w ill be overionized and recombining 
(Smit hfc Coxl l2001). In the second class of model, a se- 
ries of supernovae in a dense cloud formed a hot bub- 
ble, which burst out of the cloud into the less dense 
surroundings and underwent rap id adiabatic c ooling 
(IBreitschwerdt fc Schmutzlerlll994t IBreitschwerdt) fl996t 
IBreitschwerdt et al.l 119961: IBreitschwerdt) I2001D . In this 
case the X-ray emission is due to the delayed recombina- 
tion of the overionized ions. X-ray spectroscopy of the 
Local Bubble emission is essential for distinguishing be- 
tween the various models, as it enables us to determine 
the physical properties of the X-ray-emitting gas. 

Originally, all the soft X-ray flux was attributed to the 
Local Bubble, which made determining the Local Bubble 
X-ray spectrum relatively simple. Howev er, the discovery 
of shadows in the SXRB with ROS AT (jSnowden et all 
[19911: iBurrows fc Mendenhalll fl99lh showed that ~50% 
of the SXRB in the 1/4- keV band originated from be- 
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yond the Local Bubble, either in the Galactic halo or 
from an extragalactic background. Hence, in order to 
determine the Local Bubble spectrum, one must first dis- 
entangle the contributions of the Local Bubble and the 
background. This has been done by modeling R OSAT 
All-sky Survey data (RASS: ISnowden et al.l fl997) with 
an unabsorbed foreground component (due to the Lo- 
cal Bubble) and absorbed components for the Galac- 
tic halo and the extragalactic background, using the 
100-jnm Infrar e d Ast ronomical Satellite (IRAS) maps of 
iSchleeel et aD (fl998h as a measure of ATjj. In this way 
the Local Bubble emission has been mapped out, and 
its temperature estimated to be Tt ; r ~ 10 6 K assuming 
collis i onal ionization equilibri um (|Snowden et al.l [l998l 
[2000tlKuntz & Snowdenl[2000l) . 

It should be noted, however, that the RASS data are 
presented in just six energy bands between ~0.1 and 
~2 keV, several of which overlap. The data therefore 
have fairly poor spectral resolution. The Local Bubble 
temperature is inferred from the Rl-to-R2 band inten- 
sity ratio of the local emission component, as there is 
very little Local Bubble emission in the higher-energy 
bands: the upper limit on the ROSAT R45 intensity due 
to the Local Bubble is a few x 10 5 coun ts s" 1 arcmin -2 
()Snowden et al.lll993t iKuntz et al.lll997l ). However, this 
observational fact may be used to rule out significantly 
higher temperatures. It should also be noted that 
this inferred temperature is somewh at dependent upon 
the p lasma emission model used (e.g. IKuntz fc Snowdenl 

The large mirror collecting area and CCD detectors 
onboard XMM-Newton enable us to obtain high signal- 
to-noise spectra of the SXRB with greater spectral reso- 
lution than ROSAT. In particular, XMM-Newton enables 
us to detect line emission in the SXRB spectrum, notably 
emission from O vn and O VIII at ~0.55 and ~0.65 keV. 
If we can determine how much of this oxygen emission 
is due to the Local Bubble, this will enable us to place 
stronger constraints on Tlb- 

We use a shadowing technique to determine the spec- 
trum of the Local Bubble emission. We have analyzed 
spectra obtained from XMM-Newton pointings on and 
off an absorbing filament in the southern Galactic hemi- 
sphere. The absorbing filament appears as a shadow 
in the SXRB , as shown in Figure [TJ w hich also shows 
our X MM-Newton pointing directions. iPenprase et al.l 
(1998) have estimated the distance of the filament to be 
230 ± 30 pc. Maps of the extent of the Local B ubble, ob- 
tained from RA SS data (ISnowden et al.lll998f) and Na i 
absorption data (jSfeir et al.lll999HLallement et al.ll2003f ). 
indicate that the boundary of the Local Bubble is ^100 
pc away in this direction. Thus, the filament is between 
the Local Bubble and the Galactic halo. We fit spectral 
models simultaneously to the on- and off-filament spec- 
tra, using the difference in the absorbing column between 
the two pointing directions (iVn = 9.6 x 10 20 cm -2 versus 
1.9 x 10 20 cm -2 ; see §3. 1() to deduce the contributions of 
the foreground (unabsorbed) and background (absorbed) 
model components to the observed spectra. 

Our observations and the data reduction are described 
in Sj2j The spectral models used to fit to the data are 
described in Sj3j and the fit results are presented in §4j 
We discuss our results in <|5] In particular, we compare 
our results with the results of RASS analyses in §5.1[ and 




Fig. 1.— ROSAT All-Sky Survey R1+R2 image centered on 
£ = 279°, b = —46°, s howing the absorbing filament used for our 
observations (data from Snowdcn et al. 1997). The circles show our 
on-filament (upper) and off-filament (lower) pointing directions. 
The smaller circles (radius = 14') show the approximate areas 
from which our XMM-Newton spectra were extracted, while the 
larger circles (radius = 30') show the areas from which our ROSAT 
spectra were extracted. 

we compare our results' prediction of the Galactic halo 
O vi intensity with the observed value in £15.21 In §5.31 
we compare our results with those of other XMM-Newton 
and Chandra observations of the Local Bubble. In i j5.4l 
we discuss our choice of plasma emission code used in 
the analysis, and consider the effect this has on our re- 
sults. Finally in this section we discuss non-equilibrium 
models of the Local Bubble: an underionized model in 
£15.51 and a recombining model in ij5.6l We conclude with 
a summary in §j)J Throughout this paper we quote la 
errors. 

2. OBSERVATIONS AND DATA REDUCTION 

The on- and off-filament XMM-Newton observations 
were both carried out on 2002 May 3. The details of the 
observations are presented in Table [TJ 

Since our current understanding of the particle back- 
ground of the XMM-Newton PN is relatively poor, and 
the characterization of the background of the XMM- 
Newton MOS cameras is fairly well refined, we have re- 
stricted our analysis to the data obtained by the two 
MOS cameras. The data were reduced as follows. We 
constructed the light curve in the 2.5-8.5 keV band for 
the entire field of view. We fitted a Gaussian to a his- 
togram of the count rates, and set the "quiescent level" 
to the mean of that Gaussian. We removed from further 
analysis all time periods during which the count rate was 
>3cr above the quiescent level; the higher count rate in 
those time periods is due to either strong soft proton 
contamination or an enhanced particle background. Fil- 
tering the data using a lower-energy band (0.4-2.0 keV) 
produces results no different from those obtained using 
the above energy band. 

Sources were detected in both the 0.3-2.0 keV and 2.0- 
10.0 keV bands. Sources with a maximum likelihood de- 
tection value greater than 40 (corresponding to ^10~ 13 
erg cm -2 s -1 ) were removed. The region removed for 
each source was a circle whose radius contained 80% of 
the total flux of a point source at the source's distance 
from the optical axis; this radius was typically 24-29 
arcseconds. The few remaining faint point sources are 
likely to be background AGN with a power-law spec- 
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TABLE 1 

Details of the XMM-NEWTON observations 



Observation I b Nominal exposure Usable exposure Aoo a Nn h 

Observation ID (deg) (deg) (ks) (ks) (MJy sr _1 ) (10 20 cm" 3 ) 

On filament 0084960201 278.67 -45.32 12.8 11.9 7.10 9.6 

Off filament 0084960101 278.73 -47.09 27.8 4.4 1.22 1.9 



a 100-/im intensity from the all-sky IRAS maps oflSchlegel et al l (fl998h . b Calculated from I 100 using the conversion 
relation for the southern Galactic hemisphere in Snowdcn etTalT " ll20001 ) . 



trum which is unlikely to confuse our analysis of thermal 
spectra. The actual spectrum of the diffuse emission was 
extracted from a region with a radius of 14' approxi- 
mately centered on the optical axis after the removal of 
the point sources. 

We constructed the spectrum of the "quiescent parti- 
cle background" from the " unexposed corner" da ta and 
filter- wheel-closed data (see ISnowden et al. 2004). For 
each of our two observations, and for each MOS camera, 
the background spectrum is modeled using a database 
of filter-wheel-closed data, scaled by data from the un- 
exposed corners of the CCDs of that particular cam- 
era. This scaling is energy dependent, and is based upon 
the hardness and intensity of the unexposed corner data 
(which varies with time). The background spectrum is 
interpolated over the 1.2-1.9 keV interval before being 
subtracted from the observed spectrum. This region of 
the spectrum contains two bright instrumental lines due 
to aluminum and silicon (at ~1.48 and ~1.74 keV, re- 
spectively). In most of the spectral fits described below, 
this region of the spectrum was simply excluded. How- 
ever, leaving the instrumental lines in the observed spec- 
trum and fitting them with Gaussians during the analysis 
does not significantly affect our results. 

The strength of the residual soft proton flares is not 
known a priori, but the shape is reasonably well modeled 
by a broken power-law with a break energy of 3.2 keV, 
where the spectrum is convolved with the redistribution 
matrix but not scaled by the response function. The 
contribution of the residual soft proton flares is fitted 
during the analysis. 

2.1. Solar Wind Charge Exchange 

The solar wind was very steady during both of these 
observations. The solar proton flux, measured with the 
Advanced Composition Explorer (ACE), was ~1.8 x 10 8 
cm~ 2 s , slightly below the mean, and the proton speed 
was 420-440 km s , slightly above the mean. The 
+7 /0+ 6 and +8 /0+ 7 ratios had typical values. Both 
observations were taken at a solar angle of ^80° , avoiding 
the highest density portions of the magnetosheath. Thus, 
any solar wind charge exchange (SWCX) contamination 
will be relatively low and, more to the point, similar for 
the two observations. 

3. SPECTRAL MODELING 

3.1. Spectral Model Description 

The basic model we used to fit to our XMM-Newton 
spe ctra is based upon that use d by ISnowden et al.l (2000) 
and lKuntz &; Snowdenl (|200Q[ ) in their analyses of RASS 
data (which is i tself a development of the model used by 
ISnowden et "all 1 19981) . Thus, we use a thermal plasma 
model in collisional ionization equilibrium (CIE) for the 



Local Bubble emission. For the Galactic halo emission 
we use two thermal plasma components (2T), and for the 
extragalactic background (due to unresolved AGN) we 
use a power-law. In this model, the Local Bubble com- 
ponent is unabsorbed, while the non-local components 
(halo and extragalactic) are all subject to absorption. 

We car ried out our spectral fitting with XSPEC 1 
vll.3.2p (A rnaudl Il996f) . Our primary analysis was 
done using the A strophysical Plasm a Emission Code 
(APEC 2 ) vl.3.1 (|Smith et all l200lh for the thermal 
plasma components. For comparison, w e also ana- 
lyzed the spectra using the MeKaL model (Me we~et al.l 
Il995t ). as discussed in §5.41 For the absorption we 
used the phabs model, which uses cros s-sections from 
iBalucihska-Church fc McCammorJ (I1992D. except for H e, 
in which case the cross-section from lYaii et al.l (|1998l ) is 
used. Our basic XSPEC model was thus apec + phabs * 
(apec + apec + powerlaw). For chemical abundances 
we use d the interstellar abundance table in I Wilms et al.l 
(2000) 3 . These abundances were used both by the ther- 
mal plasma components and by the absorption model. 

The normalization and the photon index of the power- 
law used to model the extragalactic background were 
frozen at 10.5(E /keV) ~ 1A6 photons cm 2 s 1 sr 1 
keV -1 (|Chen et al.lll997t from their model A fit to ASCA 
an d ROSAT d a ta). N ote that these values were obtained 
bv lChen et afl (jl997n after removing point sources down 
to ^5 x 10~ 14 erg cm -2 s _1 , which is roughly equal to the 
limit to which we have removed sources (see @. This 
means that their result should be applicable to our anal- 
ysis. Furthermore, the exact values used for the extra- 
galactic power-law have little impact on the fitting of the 
thermal emission. 

We simultaneously analyzed the on- and off-filament 
spectra. In the fits the temperatures and normaliza- 
tions of all three apec components were free to vary, but 
were constrained to be the same for all spectra. The 
only difference in the model as applied to the different 
spectra was that the on- and off-filament spectra had 
different values of iVn for the phabs model. To deter- 
mine Nh for our two pointing directions, we obtained 
the 100-/im intensi ties for these d i rectio ns from the all- 
sky IRAS maps of lSchlegel et al.l (|1998f ) and converted 
them to ./Vh using the conversion relation for the south- 
ern Galactic hemisphere given in ISnowden et all <|2000h . 
The resulting on- and off-filament column densities are 
9.6 x 10 20 and 1.9 x 10 20 cm" 2 , respectively (see Table [TJ. 
Note that the on-filament column density is consistent 
with that derived from the color excess of the filament 

1 http: //xspec .gsfc .nasa.gov/docs/xanadu/xspec/ 

2 http: //exc .harvard . edu/atomdb/sources_apec .html 

3 Implemented using the XSPEC command abund wilm. 
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E(B - V) = 0.17 ± 0.05 (|Penprase et al.lll998l ). which 
yields -/Vh = (8.4 ± 2.5) x 10 20 cm" 2 when scale d using 
the conversion relation in lDiplas fc Savage] (|1994[ ). Also, 
our results are not very sensitive to the values of Nh 
used. 

It should be noted that our absorbing columns do 
not take into account any contribution from the warm 
ionized gas in the Reynolds layer, above the disk of 
the Galaxy. The average colu mn density of this gas 
is 7 x 10 19 / sin |6| H II cm" 2 ^Reynolds! Il99lh . which 
gives 9.6-9.8 x 10 19 H n cm -2 for our observing di- 
rections. For 1/4-keV X-rays, the effective absorption 
cross-section per hydrogen nucleus of the ionize d gas is 
62% of that for neutral gas (jSnowden et al.|[l994f ). which 
means that at low energies the Reynolds layer would ef- 
fectively contribute an extra 6 x 10 19 cm -2 to our ab- 
sorbing columns. However, we find that adding this ex- 
tra contribution to our absorbing columns affects only 
the cooler halo component, and this has no effect on our 
conclusions. Furthermore, observations of the Galactic 
absorption towards extragalactic X-ray sources are gen- 
erall y best fit without the Reyn olds layer contribution 
(e.g. lArabadjis fc Bregman1ll999f l. We therefore proceed 
using the JiMS'-derived column densities quoted above. 

As stated in Sj2] we also included a broken power-law 
component in our fit to model the contribution of resid- 
ual soft proton flares not removed in the data reduction. 
This soft proton contamination evidenced itself in certain 
XMM-Newton spectra as excess emission above that ex- 
pected from the extragalactic background power-law at 
energies > 2 keV. The broken power-law parameters 
were the same for the MOS1 and MOS2 spectra for a 
given pointing, but were allowed to differ between the 
two pointings. 

We also experimented with variants of our "standard" 
model. One variation used a non-equilibrium ionization 
(NEI) model for the Local Bubble component, i.e. we re- 
placed the first apec in our model with the XSPEC nei 
model. In this model, the emitting plasma is assumed 
to have been rapidly heated to some temperature T, but 
the ionization balance does not yet reflect this new tem- 
perature (the ions are underionized). Such a model may 
be characterized by an ionization parameter, r = n e t, 
where n c is the electron density and t is the time since 
the heating. Collisional i onization eq uilibrium is reached 
when r > 10 cm s (jMasail [19941. The nei model 
is configured to use the Astrophysical Plasma Emission 
Database (APED) to calculate the line spectrum 4 . 

Another variation replaced the 2T halo model with a 
model that used a power-law differential emission mea- 
sure (DEM) of the form 

d(E.M.[T]) _ Kjf^Y if T min < T < T max , 
d(logT) | otherwise. 

Here the exponent a and the high-temperature cut-off 
Xinax are free parameters. The low-temperature cut-off is 
frozen at T m i n = 10 5 K (the results are not very sensitive 
to the value of T m i n , as plasma at this low a temperature 
does not significantly contribute to the observed X-ray 
emission). This model is based upon the XSPEC cemekl 
model, though we modified it to use APEC instead of 

4 Implemented using the XSPEC command set neivers 2.0. 



MeKaL. The advantage of this model is that it should 
give a more accurate picture of the temperature distribu- 
tion of the gas in the Galactic halo. On the other hand, 
the 2T halo model enables an easier comparison with the 
earlier ROSAT results. 

A final variation of our "standard" model was to in- 
vestigate the effect of varying the abundances of various 
elements. 

We simultaneously analyzed the MOS1 and MOS2 
spectra obtained from the on- and off-filament pointings. 
We used the data between 0.45 and 5 keV, except for 
the region between 1.2 and 1.9 keV, where there are two 
bright instrumental lines (see ^2|. An alternative to ex- 
cluding these data is to fit these two lines by adding two 
Gaussians to the model; however, as stated in $2j there 
is no significant difference between the results obtained 
in these two ways. 

In order to better constrain the models at softer en- 
ergies, we also included two ROSAT spectra in the fits. 
These were extracted from ROSAT All-Sky Survey data 
(jSnowden et al.lll997t ) using the HEASARC X-ray Back- 
ground Tool 5 v2.3. The spectra were extracted from 0.5° 
radius circles centered on the two XMM-Newton pointing 
directions, and are normalized to one square arcminute. 
While larger circles would have reduced the errors on 
the ROSAT data, they would result in the on-filament 
ROSAT spectrum being contaminated by off-filament 
emission, and vice versa (see Fig.Q"]). 

We accounted for the difference between the XMM- 
Newton effective field of view and the solid angle used in 
the ROSAT extraction by multiplying the model applied 
to the XMM-Newton data by the XMM-Newton field of 
view (~580 arcmin 2 ), and normalizing the ROSAT spec- 
tra to 1 arcmin 2 , as mentioned above. 

After our fitting was complete, we converted the 
XSPEC model normalizations to emission measures 
(J n 2 d^) assuming tt-hc/^h =0.1, and neglecting the con- 
tribution of metals to the electron density n c . 

3.2. A Note on the Abundance Table Used 

As stated in th e previous section, we used the 
IWilms et al.l (|2000f ) interstellar abundances, which dif- 
fer f rom those in widely used solar abundance tables 
(e.g. lAnders fc Grevessdll989t iGrevesse fe Sauvallll998l ) 
for many astrophysical ly abundant elements. For exam- 
ple, IWHrnTiEaD (2000)' give an interstellar oxygen abun- 
dance ho/?ih = 4.90 x 10~ 4 , compared with so lar abun- 
dances iof_8_ 1 51_>!irrljA^ or 6.76 x 
10~ 4 (G revesse fc Sauvallll998[ ). However, more recent 
measur ements of the solar photo s pheric oxygen abun- 
dance (jAllende Prieto et all 120011 lAsplund et ai1l2004l) 
are in excellent agreement with the I Wilms et al.1 (2000) 
value. Indeed, recent measurements of other metals' 
photo spheric abundance s (e.g. C, N, Fe; lAsplund et al.1 
2005h are lower than thelAnders fc Grevessd (j!989h and 
Grevesse fc Sauval (I1998D values, a nd are in better agree- 
ment with the lWilms et ah! (|2000[) va l ues. T herefore, in 
this paper we take the IWilms et alJ (|2000f ) interstellar 
abundances to be synonymous with solar abundances. 

4. RESULTS 

4.1. Spectral Fit Results 

5 http: //heasarc .gsf c .nasa. gov/cgi-bin/Tools/xraybg/xraybg .pi 
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The results of fitting the above-described "standard" 
(2T halo) model are presented in Table [2] Also shown 
in this table are the results of using the XSPEC non- 
equilibrium nei model for the Local Bubble, and of using 
a power-law differential emission measure (DEM) for the 
halo emission (see eq. [T]). When we tried varying various 
chemical abundances, we found that the deviations from 
solar abundances were either statistically insignificant or 
inconsistent (e.g. the Local Bubble oxygen abundance 
was somewhat dependent upon the halo model used) . We 
therefore do not present any non-solar-abundance results 
in Table [21 and for the remainder of this paper we just 
discuss our solar- abundance results. 

Note from Tablc[2]tliat the Local Bubble model param- 
eters do not significantly change whether one uses a 2T 
or power-law DEM model for the halo. For completeness 
we show the halo emission measure given by the power- 
law DEM model in Figure [2j but as this paper is mainly 
concerned with the Local Bubble, we defer discussion of 
the halo results to a later paper. For the remainder of 
this paper the results discussed will be those obtained 
with the 2T halo model. 

The spectra and the best-fit "standard" model are 
shown in Figure [31 Note the offset between the total 
model emission and the extragalactic component above 
-2 keV in the off-filament XMM-Newton data. This 
is entirely due to soft-proton contamination (see §3.1j) . 
Note also that in the XMM-Newton band the Local 
Bubble component makes a significant contribution only 
to the O vii emission at -0.57 keV. Most (-90%) of 
the O vni emission at —0.65 keV is due to the hotter 
halo component, with the extragalactic background con- 
tributing —8% of the emission at this energy. However, 
the Local Bubble dominates the spectrum at the low- 
est energies in the on- filament ROSAT spectrum, as the 
halo and extragalactic components are very strongly ab- 
sorbed. The absorbing cross-section per hydrogen atom 
is 8.96 x 10~ 21 cm 2 at 0.2 keV, calcul ated using the 
iBatucinska-Church fc McCammonl (|1992f ) cross-sections 
( excep t for He; lYan et al.lll998h with the IWilms et all 
(2000) interstellar abundances. This gives an on-filament 
optical depth of 8.6 at this energy. 

While there are some features that are poorly fit in one 
of the spectra (e.g. the model underestimates the O vn 
emission in the off-filament MOSI spectrum), when one 
considers all the spectra together there are no features 
that are systematically poorly fit. Overall, the model 
gives a good fit to the data, with xl = 0-99 for 439 
degrees of freedom. 

We tested whether or not a Local Bubble component 
is necessary by trying a model without a Local Bubble 
component. As can be seen from Table [5] such a model 
gives a poor fit to the data whether one uses a 2T model 
(xl — 1-82 for 441 degrees of freedom) or a power-law 
DEM model (xl = 1.81 for 442 degrees of freedom) for 
the halo. As can also be seen from the table, a IT halo 
model gives a poor fit to the data (xl = 1-22 for 441 
degrees of freedom). 

4.2. Local Bubble O VII and O VIII Intensities 

To measure the intensities of the O vn and O vni 
emission from the Local Bubble, we replaced the Local 
Bubble apec component in our "standard" model with 
a gaussian + gaussian + vapec model, with the oxygen 
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5 5.2 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 
log (T/K) 

Fig. 2. — Our best-fitting halo power-law differential emission 
measure (DEM) model. The two arrows mark the temperatures 
obtained from the 2T halo model. 

abundance of the vapec component frozen at zero. The 
two Gaussians model the oxygen line emission, while the 
vapec component models the continuum and the contri- 
bution of other lines. The widths of the Gaussians were 
fixed at zero, and the energies were fixed at 0.5681 and 
0.6536 keV for O vii and O VIII, respectively. For O vn 
this is the mean energy of the resonance, intercombina- 
tion, and resonance lines, weighted by the line emissivi- 
ties for a 10 6 06 K plasma. For O vni it is the weighted 
mean energy of the two Lya lines. The relevant line ener- 
gies and emissivities were obtained from the APEC code 
using the XSPEC identify command. Thawing the en- 
ergies of the Gaussians did not significantly improve the 
fit, nor did it significantly alter the measured intensities. 
The temperature and normalization of the Local Bubble 
vapec component were allowed to vary, while the param- 
eters of all the other components were held fixed at their 
best-fit values from the previous section. 

We fit this new model to our on- and off-filament 
XMM-Newton and ROSAT spectra, as before. We mea- 
sure an O VII intensity of 3.4^0 4 photons cm~ 2 s _1 sr — 1 , 
and obtain a 3er upper limit on the O VIII intensity of 1.0 
photons cm~ 2 s _1 sr _1 (see Table [3]). We use these re- 
sults in £15.61 

We can check these values using line emissivity data 
from the ATOMDB database 6 . The temperature and 
emission measure of our "standard" model Local Bub- 
ble component yield O VII and O VIII intensities of 2.9 
and 0.017 photons cm -2 s" 1 sr _1 , respectively, includ- 
ing the contribution of dielectronic recombination satel- 
lite lines. These values are consistent with the above- 
measured values. Although the O Vii intensity measured 
from the Gaussian fits is higher than that obtained from 
ATOMDB, this is unlikely to be due to contamination 
by other oxygen lines, as there are no other bright oxy- 
gen lines in that vicinity. Instead, it is possible that the 
Local Bubble apec model is slightly underestimating the 
contribution of the Local Bubble to the observed O vn 
emission. This is because this component is not being 
constrained solely by the oxygen emission, but also by 
emission at other energies (e.g. the ROSAT R12 intensity 
and the Rl /R2 ratio) . The O VIII emission is more likely 
to be contaminated, for example by O vn n = 3 — > 1 
emission at 0.6656 keV. However, this just means that 
the above upper limit may be overly conservative; it does 

6 http: //cxc .harvard . edu/atomdb/download. html 



6 



HENLEY, SHELTON, AND KUNTZ 



TABLE 2 
Spectral fit results 



2T halo ! 



odcl 







Local Bubble 




Halo (cool) 


Halo 


(hot) 




Model 


logT 
(K) 


E.M. a 

(cm -6 pc) 


(10 10 cm" 3 s) 


logT 
(K) 


E.M." 

(cm -6 pc) 


logT 
(K) 


E.M. 11 
(cm -6 pc) 


X 2 /c!of 


"Standard" (CIE Local Bubble) 
NEI Local Bubble 
No Local Bubble 
One halo component 
Hot Local Bubblc c 


6-06±°° 2 

a 91 +0.08 

D - zl -o.io 

fi nc + O.Ol 
D - UO -0.02 

6.21 (frozen) 


0.018 
0.011 

0.023 
0.017 


3 4+ 4 " 7 


o.yo_ 03 

K nn + 0.02 
°- MU -0.03 

5.92 ± 0.01 
5 29+ 13 


0.17 
0.22 
0.37 

160 


6.43 ± 0.02 

6.44 ± 0.02 
6.44 ± 0.02 
6.37 ± 0.01 
6.49 ± 0.02 


0.011 
0.010 
0.012 
0.018 
0.0073 


435.86/439 
434.56/438 
801.12/441 
539.57/441 
444.04/440 






Local Bubble 






DEM halo modcl d 






Model 


logT 
(K) 


E.M. a 
(cm -6 pc) 


T b 

(10 10 cm" 3 s) 




logT max 
(K) 


a 




X 2 /c!of 


CIE Local Bubble 
No Local Bubble 


6 02+ u ul 


0.018 






6.70 ± 0.07 

>6.80 


-2 01+ 1114 

9 7 q + 0.08 
Z - ,,;> -0.07 




438.62/440 
801.64/442 



Emission measure E.M. — J n^dl. Ionization parameter r — n c t. c See j5T3] d See iTTTI for details of model. 



TABLE 3 

Local Bubble O vii and O viii intensities 



TABLE 4 

Derived parameters of the Local Bubble 



Ion 


Energy 
(keV) 


Intensity 
(photons cm" u s - sr — ) 


O VII 
O VIII 


0.5681 
0.6536 


3 4+°- 6 

< 1.0" 



a 3cr upper limit. 

n ot adversely affect o ur lat er analysis and conclusions. 

iMcCammon et all (|2002f) measured the intensities of 
O vii and O VIII in the soft X-ray background using the 
X-ray Quantum Calorimeter (XQC), which was flown on 
a sounding rocket. They obtained intensities of 4.8 ± 0.8 
and 1.6±0.4 photons cm" 2 s" 1 sr" 1 for O vn and O viii, 
respectively. These are total intensities for the soft X-ray 
background averaged over a large area of sky (~1 sr), 
whereas the intensities in Table [3] are just for the Local 
Bubble. Our results ar e therefore consistent with the 
IMcCammon et all (|2002t) results. 

4.3. Derived Parameters of the Local Bubble 

If we assume some spatial extent L for the Local Bub- 
ble in our pointing direction, and assume that the Local 
Bubble plasma is uniform along the line of sight, we can 
convert the emission measure found above to a density. 
With the measured plasma temperature, this will give 
us the thermal pressure of the plasma. If we make the 
further simplifying assumption that the Local Bubble is 
a sphere of radius L, we can estimate the thermal energy 
content and the cooling time of the Local Bubble. 

The Local Bubble parameters derived from the "stan- 
dard" model parameters in Table [2] are shown in Table [4j 

5. DISCUSSION 

In this section, we first compare our results with the 
results of analyses of ROSAT All-Sky Survey (RASS) 
data in c]5.H In c]5.2l we discuss our results in terms of 
the Galactic halo O vi emission (R. L. Shelton et al., 
in preparation). Wc compare our results with those of 
other XMM-Newton and Chandra shadowing observa- 
tions of the Local Bubble in c]5.3l We find a discrepancy 
between our Local Bubble temperature and that mea- 
sured by other XMM-Newton observations, and discuss 



Parameter 



Value 



Electron density n e (cm - ' 
Number density n (cm -3 ) 
Pressure p/k (cm -3 K) 
Thermal energy Eth (erg) 
Cooling timc a r coo i (yr) 
Sound crossing time t C ross (y r ) 



0.013(T/100 pc)~ 1/2 
0.026(L/100 pc)~ 1/2 
2.9 x 10 4 (L/100 pc)~ 1/2 
7.4 x 10 50 (L/100 pc) 5/2 
1.7 x 10 7 (L/100 pc) 1/2 
1.2 x 10 6 (L/100 pc) 



Note. — Calculated from the "standard" model parame- 
ters for the Local Bubble in Table [2] assuming the radius of 
the Local Bubble is L. 

a Calculated using the Raymond-Smith cool ing function 
([Raymond fc Smithlll977t IRaymondll 199111 with lWilms et al.l 
320001) abundances: AflO 6 06 K) = 7.5 X 10" 23 erg cm 3 s _1 . 



whether or not this can be attributed to our choice of 
plasma emission code in ij5.4l Finally we discuss non- 
equilibrium models of the Local Bubble: in c]5.5l we dis- 
cuss the results of using an underionized (ionizing) model 
of the Local Bubble (i.e. the XSPEC nei model), and in 
c]5.6l we discuss our results in terms of an overionized 
(recombining) model of the Local Bubble. 

5.1. Comparison with ROSAT Results 

In Tabic [5] we compare our measured temperatures 
with those measured in various studies of RASS data. 
As can be seen, there is excellent agreement between our 
values and the ROSA T-determined values. This is not 
surprising, as we use RASS data to constrain our spec- 
tral models at low energies. However, our Local Bubble 
emission measure (J n^dl — 0.018 cm" 6 pc) is 3-10 times 
larger than that derived from the RO SAT data (0.0018- 
0.0058 cm" 6 pc: ISnowden et al.lll998[) . 

Thi s discrepa ncy is due to the fac t thatlSnowden et alJ 
(1998) use the iRavmond fc Smith! (Il977f) plasma emis^ 
sion c ode, whereas we use APEC. Also, ISnowden et al.l 
(1998) assume a hi gher metallicity in their study: Z = 
0.017 (|Allenl 119731 ) against Z = 0.012 (jWilms et al.l 
120001 ). For a T = 10 6 06 K plasma, a Raymond & Smith 
model with Z = 0.017 predicts ~3 times as much flux in 
the 0.1-0.5 keV band as an APEC model with Z = 0.012. 
Hence, for a given amount of Local Bubble emission, our 
model will give an emission measure ~3 times larger than 
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Fig. 3. — Our observed on-filament (left-hand column) and off-filament (right-hand column) spectra, with our best-fit "standard" model 
(note the different energy ranges on the XMM-Newton and ROSAT plots). The individual components that comprise our model are also 
illustrated. Note in the on-filament XMM-Newton results that the Local Bubble (blue) and cool halo (green) components overlap. The gap 
in the XMM-Newton data is where two bright instrumental fluorescence lines have been excluded. The large difference between the XMM- 
Newton and ROSAT count-rates is because the ROSAT data have been normalized to one square arcminute, whereas the XMM-Newton 
data are integrated over the illuminated XMM-Newton field of view (~580 arcmin 2 ). 



iSnowden et alj 's (|1998f ) model, which is what we find. 

5.2. The Galactic Halo O VI Emission 

As already stated, the main focus of this paper is the 
Local Bubble emission, with discussion of the Galactic 
halo emission being deferred to a later paper. However, 
we note here that our fit results may also be used to make 
predictions of the halo O vi intensity, which may then be 
compared with that measured from an off-filament Far 
Ultraviolet Spectroscopic Explorer (FUSE) observation. 
This provides a useful check on our fit results. 

Assuming all the halo O vi emission originates from be- 
yond the absorbing material in that direction (which has 
a transmissivity to O vi photons of 58%), the intrinsic 



intensity of the doublet is 8070 



+980 
-1140 



photons cm s 



sr _1 (R. L. Shelton et al., in preparation). In compar- 
ison, the best-fit temperature and emission measure of 
the cooler halo component from our "standard" model 
predicts an intrinsic O vi doublet intensity of ^3000 
photons cm~ 2 s -1 sr _1 (using data from the ATOMDB 
database 7 ). Here we use the 2T halo model, as the halo 
differential emission measure is poorly constrained at low 
temperatures. Furthermore, we just consider the cooler 
halo component, as the contribution of the hotter com- 
ponent to the O vi emission is negligible. 

We used a Monte Carlo method to estimate the un- 
certainty on this O vi intensity prediction. We gener- 
ated 1000 random pairs of (T [halo], E.M. [halo]), where 

7 See footnote [6] 
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TABLE 5 
Comparison with ROSAT results 



Work 



log T LB 
(K) 



logTW 
(K) 



log T H al< 

(K) 



Snowden et al. f 1998*i a 
Snowden et al. (2000) 
Kuntz fc Snowden (20001 
This work d 



6.07 ± 0.05 6.02 ±0.08 



6.08° 
6 ll+ 015 

q + 0.02 

-0.04 



6.06Z 



6.00° 
6.06I 20 

, + 0.04 
-0.03 



5.93Z 



6.4° 
6 46+ 012 
6.43 ± 0.02 



a No second halo component used. No errors quoted. c Temperature 
fixed at this value. d "Standard" model results from Table [2] 
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Fig. 4. — Histogram of 1000 randomly generated values of the 
intrinsic Galactic halo O VI intensity, predicted from our measured 
values of the halo temperature and emission measure, taking into 
account the errors on these parameters. The arrow and the hor- 
izontal bar denote the intrinsic intensity inferred from FUSE, as- 
suming all the O VI emission originates from beyond the absorbing 
material in that direction (see text for details). 

E.M. [halo] = J ' n^dl is the emission measure of the 
cooler halo component. These pairs of numbers were 
drawn at random from normal distributions whose stan- 
dard deviations are given by the errors on the measured 
parameters in Table|l(E.M. [halo] = 0.17±g;^| cm" 6 pc), 
and were used to calculate 1000 values of the O vi dou- 
blet intensity. We fit a Gaussian to the distribution of 
these values, and from the mean and standard deviation 
obtain a predicted O vi intensity of 3100 ± 1000 photons 
cm -2 s _1 sr -1 (see Fig.[4]). The observed intensity is 3.3a 
larger than the predicted intensity (though the distribu- 
tion of predicted intensities is positively skewed, which 
will tend to reduce the significance of this difference). 
These results indicate that there is more O VI in the halo 
than is expected from the hot gas alone. This is prob- 
ably because O vi can also arise in the warm interfaces 
between cool clouds and the hot gas. 

5.3. Comparison with Other Shadowing Observations of 
the Local Bubble 

XMM-Newton has been used to make shadowing 
observations of the Local Bubble in the directions 
of the MBM 12 and O phiuc hus molecular clouds 
(|Frevberg fc Breitschwerdtl 120031 : iFrevberd l200l . and 
in the direction of the Bok globule Barnard 68 
(|Frevberg et al.ll2004|). Chandra ha s also been used to 
observe MBM 12 dSmith et"aT][2005h . 

The Local Bubble temperatures inferred from these 
XMM-Newton observations are consistently higher than 
the temperature we measure: log(TLB/K) = 6.2lj^;°y 
for MBM 12 and Ophi uchus (|Frevberg k, Breitschwerdtl 
l200l IFrevberd l200l and log(T LB /K) k, 6.24 for 



Barnard 68 (see Appendix), versus log(TLB/K) = 
6.061 

ao4 fr° m our data. Note, however, that for 
MBM 1 2 there is an additional log(T/K) sa 6 com- 
ponent (iFrevber d I2004D. It sho uld also be noted t hat 
iFrevberg fc Breitschwerdtl (|2003D and lFrevbergl 1)2004 ) do 
not state whether or not they use ROSAT data to con- 
strain their fits at lower energies. 

In c ontrast to these XMM-Newton results. ISmith et"ail 
(2005) found that they could not explain the O VII:0 vm 
ratio in their Chandra observation of MBM 12 with 
an equilibrium Local Bubble model with log(TLB/K) < 
6.3. However, they do note the possibility that their 
O vm emission is contaminated by emission from an- 
other source, such as solar charge exchange emission. 

These other shadowing observations were carried out 
using much thicker absorbers than our observa tions: 
4 x 10 21 cm" 2 for MBM 12 (ISmith et all [20051) u p 
to — 10 23 cm -2 for Barnard 68 (|Frevberg et all I2004D . 
against 9.6 x 10 20 cm" 2 for our filament. The optical 
depth at O vm energies is at least 2 for these other 
observations, implying a transmissivity for background 
O vm radiation of less than 14%. In an observation of 
one of these thicker absorbers, one may confidently at- 
tribute a large fraction of any observed O vm emission 
to the foreground, and thus infer a higher Local Bub- 
ble temperature. However, the higher transmissivity of 
our filament (61% at O vm energies) makes it harder 
to determine how much of the observed emission is from 
the Local Bubble, and how much is background emis- 
sion that has leaked through the filament. Our best-fit 
"standard" model attributes only ~2% of the observed 
on-filament O vm emission to the Local Bubble, com- 
pared with ~30% of the observed O vn emission. The 
fact that our best-fit model attributes so little O vm 
emission to the Local Bubble leads to our lower value of 
Tlb- However, the question remains, could more of the 
O vm emission in our observation be due to the Local 
Bubble? To put this more precisely, are our data also 
consistent with log (Tlb /K) = 6.21? 

We tested this by re-fitting our "standard" model to 
the data, but with T LB frozen at 10 6 21 K. The results 
of this are shown in Table [2j This model does give 
a good fit to the data (x 2 = 1-01 for 440 degrees of 
freedom). However, note that the temperature of the 
cooler halo component has significantly decreased, and 
its emission measure has increased by three orders of 
magnitude. Both these effects lead to a huge increase 
in the predicted intrinsic halo O VI intensity (cf. £15. 2|) 
to ^10 7 photons cm -2 s -1 sr -1 . As before, we used 
a Monte Carlo method to estimate the uncertainty on 
this prediction. In this case, there is a much larger error 
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on the emission measure of the cooler halo component 
(E.M. [halo] = lGOljo 00 cm ~ 6 P c )> resulting in a much 
larger dynamic range in the predicted intensities. By fit- 
ting a Gaussian to the distribution of the logarithms of 
the predicted values, we find the predicted intensity is 
log(/ovi/ph cm" 2 s _1 sr _1 ) = 7.5 ± 1.1. In comparison, 
the intrinsic intensity inferred from the FUSE data is 
log(Jovi/ph cm" 2 s" 1 sr" 1 ) = 3.79±g;gj| (R. L. Shelton 
et al., in preparation). While we could interpret the dis- 
crepancy between the predicted and observed intensities 
as evidence that the Galactic halo is out of equilibrium, 
given the size of the discrepancy (i.e. possibly up to a 
few orders of magnitude), we instead interpret it as ev- 
idence that log(TLB/K) = 6.21 is inconsistent with our 
XMM-Newton spectra and the FUSE results when taken 
together. 

It should be noted that two of the clouds discussed 
above (Ophiuchus and Barnard 68) lie beyond the Lo- 
cal Bubble in the direction of Loop I, near the Galactic 
plane. It is therefore possible that the foreground emis- 
sion (which in the above discussion we have attributed 
to the Local Bubble) is being con t amin ated by emission 
from Loop I. However, iFreyber j (|2004f ) states that the 
weakness of the Fe-L line emission (attributed to Loop I) 
in the on-cloud Ophiuchus spectrum indicates that the 
contamination is small. In contrast to this, MBM 12 does 
not lie towards any obvious source of contamination, and 
may even lie within the Local Bubble, implying that the 
temperature inferred from these observations is indeed 
that of the Local Bubble. 

Our pointing direction is > 80° away from the direc- 
tions of the other XMM-Newton observations. It is there- 
fore not implausible that our value of Tlb is different 
from those measured from these observations, as we are 
observing a different part of the Local Bubble. These re- 
sults therefore suggest that the Local Bubble is thermally 
anisotropic. 

5.4. Choice of Plasma Emission Code 

In this section we consider what effect, if any, the choice 
of plasma emission code used has on the measured Local 
Bubble temperature. 

The higher foreground temperatures measured from 
the XMM-Newton spectra of MBM 12, the Ophiuchus 
molecular cloud, and Barnard 68 should be quite robust, 
regardless of the plasma emission code used. This is be- 
cause O VIII emission is observed towards these clouds. 
Due to the large optical depths of the clouds, a large 
fraction of this emission may be attributed to the fore- 
ground, implying a hi gher temperature than we found 
towards our filament. iFreyberg fc Breitschwerdn (2003) 
do not give details of the models used in their analysis 
of the XMM-Newton observations of MBM 12 and Ophi- 
uchus, but they do note that log(T LB /K) = 6.21 
"or higher, depending on the a ctual model." Al so, we 
infer log(T LB /K) = 6.24 from IFreyberg et all 's (|2004t ) 
Barnard 68 data whether we use APEC or MeKaL (see 
Appendix). However, since our filament is less optically 
thick than these other clouds, there is a greater ambigu- 
ity between what is Local Bubble emission and what is 
halo emission, and so we test the possibility that different 
codes would attribute different amounts of the emission 
in our spectra to the Local Bubble and the halo, leading 
to a different value of Tlb- 
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Fig. 5.— Our best fitting APEC and MeKaL Local Bubble 
models, folded through the XMM-Newton MOS1 instrumental re- 
sponse. Note that there is more Local Bubble O VIII emission at 
~0.65 keV in the MeKaL model than in the APEC model. 

We tested this by repeating our fits using the MeKaL 
code instead of APEC for the thermal emission compo- 
nents. The results are compared with our APEC results 
in Tabled Note in particular that MeKaL gives a higher 
Local Bubble temperature than APEC: log(TLB/K) = 
6.171^7 (MeKaL) versus 6.061^4 (APEC). It should 
be emphasized that the difference between the APEC 
and MeKaL results for our data is not because the 
two codes give different temperatures for the same Lo- 
cal Bubble spectrum. Instead it is because MeKaL at- 
tributes more of the O vin emission to the Local Bubble 
(see Fig. [5J, and correspondingly less to the halo. Hence, 
the inferred Local Bubble spectrum is different between 
the two codes, and thus so too is Tlb- This discrepancy 
is most likely to be due to uncertain modeling of the 
lines from L-shell ions of Ne, Mg, and Si, which domi- 
nate the emission at the lowest ROSAT energies. The 
ATOMDB vl.3.1 release notes 8 contain a caveat that 
there are very few data on lines from these ions (other 
than Li-like ions). Differences between the codes in this 
energy regime would affect the fits to the ROSAT data, 
which would then affect the fitting to the higher-energy 
XMM-Newton data. 

To test whether or not this does cause the discrep- 
ancy between the APEC and MeKaL results, we re-fit 
our "standard" model just to the XMM-Newton spectra, 
without the ROSAT data. These results are also shown 
in Table O Note that these temperatures are unphysi- 
cal, as they are inconsistent with the low-energy ROSAT 
data. However, the agreement between the codes' results 
is much better than before, suggesting that at least part 
of the discrepancy is due to uncertain modeling of lines 
in the lowest-energy ROSAT bins. 

Despite this discrepancy, we reiterate that the mea- 
surement of log(T LB /K) > 6.21 from the other XMM- 
Newton observations should not be strongly code- 
dependent, and note that we obtain a lower tempera- 
ture than this whether we use APEC or MeKaL. This 
therefore implies that the conclusion of the previous sec- 
tion, namely that the Local Bubble appears to be ther- 
mally anisotropic, is not an artefact of our choice of 
plasma emission code. However, we should reiterate the 
possibility that the foreground emission in some of the 
other XMM-Newton observations may be contaminated 

8 http: / /exc .harvard . edu/atomdb/ issues_improvements .html 
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by non-Local Bubble emission. 

5.5. Ionizing Model of the Local Bubble 

In the preceding sections we have just discussed col- 
lisional ionization equilibrium (CIE) plasma models. In 
this section we discuss the results of using the XSPEC 
nei model to model the Local Bubble emission (see Ta- 
ble[5]). This is a model in which the ions are underionized, 
i.e. the plasma has been rapidly heated, but the ioniza- 
tion balance does not yet reflect the new temperature: 
the ions are still in the process of ionizing. As stated in 
£13.11 such a model may be characterized by an ionization 
parameter, r = n e t, where n e is the electron density and 
t is the time since the heating. Eq uilibrium is reached 
when r > 10 12 cm" 3 s (|Masaiill994l) . 

In certain situations one may use the F-test to de- 
termine whether or not adding an extra model pa- 
rameter leads to a signifi cant improvement in \ 2 
(|Bevington fe R obinson 120031 ). Unfortunately we can- 
not use the F -test here to test whether or not the non- 
equilibrium model is an improvement on the equilibrium 
model. This is because the additional parameter (r) is 
on the boundary of the set of possible parameter values 
in the simpler (null) model (i.e. r = oo in the equilib- 
rium model). In such a case, one cannot use the -F-test 
(|Protassov et al.ll2002t) . 

However, we can note that the difference in \ 2 between 
the two models is very small, suggesting that thawing r 
does not lead to a significant improvement in \ 2 (even if 
we cannot use the F-test to demonstrate this for certain) . 
On the other hand, it should be noted that the value of r 
we measure ([3. 4^4] x 10 10 cm -3 s) is significantly less 
than w hat is expect ed in the equilibrium case (r > 10 12 
cm" 3 s; Masai 1994), which means that our data cannot 
rule out an ionizing non-equilibrium model. 

In order to distinguish between equilibrium and non- 
equilibrium models we would require more spectral in- 
formation. With our data, the only Local Bubble "line" 
(actually a blend of lines) we can clearly see is the O vn 
feature at ~0.57 keV. There are no other obvious Lo- 
cal Bubble lines in the XMM-Newton bandpass (the ob- 
served O Viii feature in the on-filament spectrum is at- 
tributed to leak through from the Galactic halo), and 
we only have two ROSAT spectral bins at energies be- 
low the XMM-Newton bandpass. In contrast to this, 
one needs to compare the flux ratios of lines from sev- 
eral different ions in order to determine whether or not 
non-equilibrium effects arc important. 

In principle, we could use data from other wavebands 
to help distinguish between the two classes of model. 
FUSE has observed the Local Bubble in the same di- 
rections as our XMM- Newton observa tions. Using the 
on-filament observation, ISheltonl (|2003[ ) has placed a 2a 
upper limit on the Local Bubble O vi AA1032, 1038 dou- 
blet intensity of 800 photons cm" 2 s" 1 sr" 1 (including 
statistical and systematic uncertainties). In compari- 
son, the predicted O vi intensities of our best-fit models 
(determined using XSPEC) are 180 photons cm" 2 s" 1 
sr" 1 (equilibrium) and 63 photons cm" 2 s" 1 sr" 1 (non- 
equilibrium) , both of which are consistent with the ob- 
served upper limits. 

We therefore conclude that our data are unable to dis- 
tinguish between an equilibrium plasma model and a 
non-equilibrium model in which the ions are underion- 



ized and are in the process of ionizing. 

5.6. Recombining Model of the Local Bubble 

In this section we consider an alternative non- 
equilibrium Local Bubble model to the previous section, 
namely one in which the plasma is overionized and the 
ions are in the process of recombining. 

Such a model for the Local Bubble was origi- 
nally pr oposed by Breitschwcrdt & Schmut zlerl ((19941 ; 
see also IBreitschwe rdt 1996; B reitschwerdt et "all 119961 : 
iBreitschwerdtj 120011 ) . In this model, the Local Bubble 
is assumed to be the relic of an old superbubble formed 
by ~10 successive supernovae in a dense (~10 4 cm" 3 ) 
molecular cloud. A few million years ago, the bubble 
burst out of the cloud and underwent rapid adiabatic 
expansion into the surrounding medium, assumed to be 
~10-100 times less dense than the cloud. During this 
process the bubble adiabatically cooled to < 10 5 K (i.e. 
much less than the temperature inferred from the X-rays 
assuming collisional ionization equilibrium). However, 
owing to the rapidity of the cooling, the ions remained in 
high ionization states. The observed X-rays are then sup- 
posed to be due to the delayed recombinations of these 
overionized ions, rather than being from a hot (~10 6 K) 
plasma in equilibrium. 

A recombining plasma model has a number of appeal- 
ing features. Firstly, there is an order of magnitude 
difference between the pressure of the Local Cloud in 
which the solar system resides {Phc/k ~ 2000 cm" 3 K; 
lLallementlll998l ) and the pressure of the surrounding Lo- 
cal Bubble inferred from the X-ray data assu ming CIE 
(PhB/k ~ 15,000 cm" 3 K [Snowden et aflll998l ] to 29,000 
cm" 3 K [Table H]). It is difficult to see how such a large 
pressure difference can be maintained. The lower Local 
Bubbl e temperature in the Breitsch werdt fc S chmutzler 
(|1994f ) model greatly reduces the Local Bubble pres- 
sure, and thus eliminates this problem. Secondly, the 
electron density inferred from the dispersion measure of 
the pulsar PSR 0950 + 08 (d « 130 pc) is 0.023 cm" 3 
(|Revnoldslll9"90l) . If this density is representative of the 
Local Bubble, then a ~10 6 K plasma in equilibrium 
would produce too much X-radiation (cf. the electron 
density inferr ed from the X-ray da ta assuming CIE is 
~0.007 cm" 3 [Snowden et aD 119981 ] to 0.013 cm" 3 [Ta- 
bic |U ). However, a higher n f is a cceptable within the 
IBreitschwerdt fc Schmutzlerl (| 19941 ) model. Thirdly, a 



~10 6 K plasma in equilibrium should produce most of 
its emission as lines in the extreme ultraviolet (EUV). 
However, such emission is not observed. Extreme Ultra- 
violet Explorer (EUVE) spectra of the diffuse EUV back- 
ground (~ 17-80 eV) place an upper limit on the emission 
measure of a local 10 6 K plasma which is an order of mag- 
nitude l ess than the Local Bubble emission measure s in 
Table l2l ([ Jelinksv et al.lfl995t iVallerga fc Slavinlll99l . A 
recombining plasma model can explain this observation, 
as the much lower kinetic tempera ture suppresses coll i- 
sional excitation of the EUV lines ([Breits chwerdt 1996). 

Here we examine whether or not a recombining plasma 
model can simultaneously explain our observed Local 
Bubble O vn intensity and the observed upper limits on 
the O VI and O Viii intensities. We do this by first cal- 
culating the fractions of oxygen in the ionization states 
O vi to O ix, and then using these populations to pre- 
dict the line intensities. These are then compared with 
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TABLE 6 

Comparison of APEC and MeKaL results 



Code 


ROSAT 
included? 


logT LB 
(K) 


logT H alo,l 

(K) 


logT Ha lo,2 

(K) 


X 2 /dof 


APEC 
MeKaL 
APEC 
MeKaL 


Y 
Y 
X 
N 


6 06+ 02 

u - uu -0.04 

6 17+ - 013 

011 -0.07 

6 44+0-03 

"■^-0.08 
O" 3 ' -0.10 


K QQ + 0.04 
°- ao -0.03 
r fi7 +0.01 
°- DI -0.02 

6 20+ 04 
6 21 +0 13 


6.43 ± 0.02 

6.35 ± 0.02 

6 69+ ' 10 

6 7^+°- 15 
D - ,o -0.20 


435.86/439 
433.68/439 
410.21/425 
411.57/425 



the obser ved limit s . For this purpose we use 3a limits. 
Note that ISheltonl (|2003l ) does not explicitly quote a 3a 
upper limit for the O vi intensity, but we can obtain one 
from the data in her Table 2. The 3a upper limit on the 
A1032 intensity obtained from the combined day + night 
FUSE data is 710 photons cm -2 s" 1 sr _1 (statistical 
uncertainty only). As the theoretical A1032-to-A1038 in- 
tensity ratio is 2 to 1, the 3a upper limit on the doublet 
intensity is 1065 photons cm -2 s -1 sr _1 . 

In a purely recombining plasma (i.e. one in which ion- 
izations can be ignored), the fraction pi of the atoms of 
a given element in an ionization state with i electrons is 
given by 

d^?i 

— = -RiUePi + Ri-in e pi-i (2) 

where Ri is the recombination coefficient (radiative + 
dielectronic) to go from the ?th state to the (i + l)th 
state (e.g. Rq is the recombination coefficient for O ix 
— > O vni). Note that for i = there is no second term 
in equation ([2]), and it may easily be solved for po- This 
solution m ay then be used to solve equation ^fy for p\ , 
and so on. ISmith et alJ (2005) give the solutions for this 
system of coupled differential equations for the case of a 
plasma that is initially fully ionized, i.e. po = 1, Pi>o = 
at t = (their eqs. [5] and [6]). We used these solutions 
to calculate the populations of O ix to O VI, using ra- 
diative recombination coefficients from lVerner fc Ferlandl 
(119961) (the excep tion to this is O vi — > O V, for which 
iVerner fc Ferland| do not give data; in this case we use 
data from lShull fc van Steenbergjll982f). and die l ectron ic 
recombination coefficients from Mazzott a~et al.1 JI998). 

To calculate the O vn and O vni intensities, we use 
rate coefficients a R for the formation of the lines of inter- 
est due to recombinations (radiative + dielectronic) from 
the above ionizati on stage, calculated using data from 
IMewe et al.1 (|I985). We can neglect the contribution of 
collisional excitations to the lines, as the temperature 
of the cooled plasma (kT < 20 eV) is much less than 
the energy required to excite an electron to the n = 2 
level (~570 eV). The intensity lovn of O vn emission 
is therefore given by 

af f 

I ovu = — / n c n viudl 

= (tr)Gl)M < 3 > 

where of- is the rate coefficient for the production of 
O vn Ka (i.e. n = 2 — * 1) due to recombinations from 
O VIII, comprising contributions from the resonance, in- 
tercombination, and forbidden lines (the contribution 
of unresolved satellite lines is negligible in this model; 
IMewe fc Gronenschildlll981fh In equation © we have 



used 

/ no vm \ / n \ ( nu\ 
no vni = — — 

V n J \ n H/ V n e J 

where no vm/ n o = Pi an d %/n c = 1/1-2 for a fully ion- 
ized plasma with njj/n-Hc = 10- Similarly, the intensity 
-fo vni of O VIII emission is given by 

where aj^ is the rate coefficient for the production of 
O VIII Lya due to recombinations from O IX. 

For O VI, collisional excitation dominates over recom- 
binations from O vn. The intensity of the AA1032, 1038 
doub let may be calculated by rearranging equation (5) 
from lShull fc Slavinl (jl99l : 

O(T) / 1.392A 
/ OVI = _exp(-— j 

where T5 is the gas temperature T in units of 10 5 K, 
and f2(T) is the Maxwellian-averaged collision strength 
for de-excitation of the O vi doublet. We use the 
parameterizati o n for H(T) given by equation (6) in 
IShull fc Slavinl (|1994f K We calculate novi from n c us- 
ing the oxygen abundance and P3 = novi/no- ^ovi — 
P3(no/n H )K/1.2) (cf. cq. g]). 

We have calculated oxygen line intensities as functions 
of the time t since recombinations began for electron den- 
sities between 0.001 and 0.3 cm -3 and temperatures be- 
tween 10 4 and 10 7 K. Figure [6] shows an example of our 
res ults, calculated for n P . = 0.024 cm" 3 and T = 4.3 x 10 4 
K (|Breitschwerdtll200lD . For a given (n e , T), the tightest 
constraint on t is set by the measured O VII intensity. 
The predicted O vi and O vni intensities at this time 
may then be compared with the observed upper limits. 
For example, in Figure [5] the observed O vn intensity 
implies that t ~ few x 10 3 or ~3-4 x 10 5 yr. The former 
case is strongly ruled out as the predicted O vni inten- 
sity is an order of magnitude larger than the observed 
3a upper limit. In the latter case, however, both the 
O vi and the O Vin intensities are within the observed 
3a upper limits. 

For each (n e , T) we investigated, we check whether or 
not there exists a time (more precisely, a range of times) 
for which all three predicted line intensities lie within 
their observed 3a limits. The results are illustrated by 
the black contour in Figure [7] - in the region above this 
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Fig. 6. — O VI, O VII, and O vm intensities predicted by our 
recombining plasma model as a function of time (the plasma is fully 
ionized at t = 0). The curves were calculated using n c = 0.024 
cm" 3 and T = 4.3 X 10 4 K l!Breitschwerdt|[200ll) . The horizontal 
lines show our observed Local Bubble O VII intensity (with the 3<r 
uncertai nty; see T abled, and the 3<r upper limits on the intensities 
of O VI (Shelton 2003) and O vm (see Table©. 




Fig. 7. — Black contour: for values of (n e , T) below this con- 
tour there exists a time for which all three predicted line intensities 
simultaneously lie within the observed 3cr limits. Above this con- 
tour no such time exists. Red contours: contours showing the time 
t after recombinations began at which the predicted O VII inten- 
sity is closest to the observed value (from bottom to top, t = 6, 
7, and 8 X 10 5 y r). Black cross: This marks n c = 0.024 cm , 
T = 4.3 X 10 4 K HBreitschwerdtil200l]; see Fig. ©, 

contour, no time exists for which all three lines are si- 
multaneously within their observed limits. In the left 
portion of this region, the O vn intensity nevers reaches 
the observed value, while in the right portion, either the 
O vi or the O vm is too bright whenever the O vn is 
within the observed limits. We note that this is not a 
rigorous statistical test. Such a test would be difficult, if 
not impossible, as we only have three data points, two of 
which are upper limits (so, for example, the \ 2 test can- 
not be used). However, it does help to illustrate which 
regions of (n e , T)-space are and are not observationally 
acceptable. 

Also for each (n c , T) we find the time for which the pre- 
dicted O vn intensity is equal to the observed value (or 
closest to it, if it never actually equals it). This provides 
the tightest constraint on the time t since recombina- 
tions began. If there are two such values (as in Figure[5]) 
we take the later one. These results are illustrated by 
the red contours in Figure [71 As can be seen, the obser- 
vationally acceptable region of (n e , T)-space (i.e. below 
the black contour) corresponds to t < 7 x 10 5 yr, and 
most of this region (including the values of n e and T 



from iBreitschwerdtl l2001h corresponds to t < 6 x 10 5 
yr. Furthermore, the observationally acceptable region 
above the 6 x 10 5 -yr contour corresponds to T > 5.5 x 10 5 
K and n c > 0.018 cm -3 . This gives a pressure P > 20000 
cm~ 3 K, which is an order of magnitude lar ger than the 
pressure of the Local Cloud l)Lallementlll998f ). While this 
is not as strong a constraint as the observed oxygen line 
intensities, it does mean that in this region of parameter 
space the recombining model loses one of its main selling 
points. 

In summary, for most of the observationally accept- 
able region of (n c , T)-space, the observed O vn inten- 
sity implies that the time since recombinations began is 
t < 6 x 10 5 yr. We shall now show that this time is also 
an upper limit on the time since the Local Bubble burst 
out of its natal cloud, as the timescale for recombinations 
was shorter when the Local Bubble was smaller. 

If the Local Bubble is expanding and cooling adiabat- 
ically, its temperature T and volume V are related by 

TV 1 ' 1 = constant, (7) 

where 7 is the adiabatic index (7 = 5/3). Therefore, the 
temperature of the Local Bubble when it is of linear size 
L is 

T = T Q (L/L y 2 , (8) 

where To and Lq are the present-day temperature and 
size. Similarly, the electron density n e is related to the 
present-day density n c0 by 

n c = n c0 (L/L Q )- 3 . (9) 

The recombination timescale t rcc is defined by 

^rcc = r>fr r \ 1 (10) 

where R(T) is the recombination coefficient (see eq. 0). 
This timescale may be expressed as a function of L/Lq, 

t rcc _ (L/L ) 3 R(T ) 
t rec ,o i?(T [L/L ]- 2 )' 1 ' 

where t rec is the present-day recombination timescale. 

Figure^a) shows t Iec /t TeCi o for O VI, O VII, and O VIII 
as functions of L/Lq, assuming a present- day tempera- 
ture T = 4.3 x 10 4 K (|Breitschwerdtl[200lb . Figure^b) 
shows the recombinations and ionization coefficients for 
the various process that govern the evolution of the pop- 
ulation of O VIII, again as functions of L/Lq. The 
ionization coefficients were calculated using data from 
lArnaud fc Rothenflud (| 19851 ). 

From Figure Ufa) we can see that the recombination 
timescale was always shorter at earlier times, when the 
Local Bubble was smaller. This is true for all present-day 
Local Bubble temperatures we investigated. From Fig- 
ure Efb), meanwhile, we can see that when L > 0.15-Lrj, 
the recombination coefficients governing the evolution of 
the population of O VIII are larger than the ionizat i on co- 
effic ients (cf. in the iBreitschwerdt fc Schmutzlerl [1994| 
and iBreitschwerdt et alJ |1996l | models, the Local Bub- 
ble burst out of its natal cloud when L ps 0.2L )- This 
means that since the Local Bubble burst out, recombi- 
nation has dominated over ionization in the evolution of 
the O vm population. We find this is true in general for 
Tq < 10 s K. This result, combined with the fact that the 
timescale for recombinations was shorter when the Local 
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Fig. 8. — (a) Recombination timescales t Tec for the formation 
of the labeled ions as a function of Local Bubble size L, assuming 
the Local Bubble is expanding adiabatically. The recombination 
timescales and L have been normalized to their present-day val- 
ues t rec ,o and Lo, respectively, (b) Recombination and ionization 
timescales for the labeled processes governing the evolution of the 
population of O VIII. In both figures the curves were calculated 
assuming a present-day temperature To = 4.3 X 10 4 K. 

Bubble was younger, means that the population of O vni 
will be depleted more rapidly than in our original model 
in which n e and T are held fixed at their present-day 
values. 

The discussion in the previous paragraph applies only 
t° Tq < 10 5 K. For higher present-day temperatures 
(To ~ few x 10 5 K), the ionization coefficients that gov- 
ern the O VIII evolution remain larger than the recom- 
bination coefficients for some time after the Local Bub- 
ble has burst out of its natal cloud. However, even for 
Tq = 3 x 10 5 K, recombinations start to dominate the 
evolution of O vni before the Local Bubble reaches half 
its current size. Furthermore, the much lower recombi- 
nation timescale at this time more than compensates for 
the delay in the start of recombinations - by the time the 
Local Bubble reaches its current size, the O vni popula- 
tion will be lower than if the ion populations had evolved 
for the same amount of time with n Q and T held fixed at 
their present-day values. 

The result of the above discussion is that the O vn 
emission (which in this model is formed by recombina- 
tions from O vni) will fall below the observed level on a 
timescale even shorter than the ^6 x 10 5 yr given by our 
original model. Hence, as stated above, 6 x 10 5 yr is an 
upper limit on the time since the Local Bubble burst out 
of its natal cloud. If the Local Bubble were any older 
than this, the O vn emission would be below what is 
observed. 

This age is a factor of 2-7 times smaller than the age 
of the Local Bubble proposed in recombining models 
of the Local Bubble (IBreitschwerdt & Schmutzlerlll994t 
iBreitschwerdt et ail I1996T ). It also means that the Lo- 



cal Bubble would have expanded extremely rapidly after 
bursting out of its natal cloud: (100 pc)/(6 x 10 5 yr) rj 
160 km s _1 . This is much larger than the adiabatic sound 
speed of the ambient medium (c s = 15 km s _1 when 
T = 10 4 K), and so the expansion would have been highly 
supersonic, resulting in a shock at the edge of the Local 
Bubble as it expanded, with associated heating and ion- 
ization processes. Thus, the above-described model of an 
adiabatically cooled, purely recombining Local Bubble is 
not self-consistent. 

It should be noted that strictly the above argument 
against a recombining model of the Local Bubble only ap- 
plies for present-day temperatures To < 3 x 10 5 K. How- 
ever, as noted above, higher present-day temperatures 
which are observationally acceptable (see Fig. [7]) imply a 
Local Bubble pressure an order of magnitude larger than 
that of the Local Cloud, in which case the recombining 
model loses one of its major selling points. It is also 
possible that such higher temperatures would still lead 
to an implausibly young Local Bubble age. However, to 
test this would require more detailed calculations of the 
growth of the Local Bubble and the ionization evolution 
therein, which are beyond the scope of this paper. 

This is not the first time objections have been raised 
against a re combining model of t he Local Bubble. Using 
values from IBreitschwerdt! (|2001l ). and taking as a lower 
limit for the O vi fraction the value f or a 5 x 1 5 K p lasma 
in equilibrium (novi/ n o = 0-037). ISheltonl (|2003l ) pre- 
dicted an O vi intensity of ^1900 photons cm -2 s^ 1 sr _1 . 
She adds that this prediction is a lower limit, as the O vi 
fraction would tend to increase due to recombinations 
from O vn, and would not decrease again until most of 
the oxygen had recombincd to O I and O II. This pre- 
dicted intensity is more than twice the 2a upper limit 
on the Local Bubble intensity established with FUSE 
(|SheltorJ l2003f). T his, and a number of other discrep- 
ancies, led ISheltonl to conclude that recombining models 
c ould be practically el iminated. 

lOegerle et al.l (|2005f ) found that a recombining model 
of the Local Bubble also disagrees with observations from 
the point of view of O vi absorption. They found that 
such a model leads to an O VI column density much larger 
t han that which they had measured with FUSE. 

ISmith et all 0m \ considered a large set of recombin- 
ing plasma models, covering a wide range of electron 
densities, temperatures, and Local Bubble sizes. They 
found that they were unable to match simultaneously 
the O vn and O vm intensities they measured from their 
Chandra spectr um of MBM 12, the upper limit on the 
O VI intensity (jSheltonl 120031 ). the O VI column density 
(|Oegerle et al.ll2005h . and the ROSAT R12 intensity. In 
essence, the R12 intensity requires a minimum density 
of highly ionized ions; however, the densities of O VI, 
O vn, and O vm are limited by the F USE and Chandra 
data. Given this, ISmith et alJ (|2005l ) said that a static 
recombining plasma model (i.e. a model in which the 
temperature and density governing the recombinations 
is const ant) may b e confi dentl y disposed of . 

Both |Sheltoi's (|200l and (Smith etaLT s (l200l ob- 
jections to a recombining model are based upon the 
predictions of a sta tic recombining model. However, 
IBreitschwerdt! (|2001h points out that for a proper re- 
combining model of the Local Bubble, the dynami- 
cal and thermal evolution of the Local Bubble must 
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be treated together, self-consistently. While we have 
not carried out full dynamical modeling of the ion- 
ization evolution in an expanding, cooling Local Bub- 
ble, we have shown that our static model gives a 
Local Bubble age several times smaller than th at in 
the models of IBreitschwerdt fc Schmutzlerl (| 19941 ) and 
iBreitschwerdt et al.l ( 19961 ), and furthermore by consid- 
ering the recombination timescales when the Local Bub- 
ble was smaller, we have shown that the Local Bubble 
age given by our static model is an upper limit on the 
age a dynamical model would give. To put this another 
way, if the Local Bub ble had burst out of its natal cloud a 
few million years ago (IBreitschwerdt fc Schmutzlerill994t 
IBreitschwerdt et al.lll996l) . the O VII intensity would be 
far below what is observed. Our observations therefore 
add to the evidence against a recombining plasma model 
(static or dynamical) of the Local Bubble (at least with 
a present-day temperature <3x 10 5 K). 

6. SUMMARY AND CONCLUSIONS 

We have analyzed XMM-Newton spectra of the inter- 
stellar medium, obtained from pointings on and off an 
absorbing filament at high southern Galactic latitude. 
We have fit various models simultaneously to both sets of 
spectra, and used the difference in the absorbing column 
in the two pointing directions to constrain the spectrum 
of the Local Bubble. 

Our main findings are as follows: 

1. In the examined direction, The Local Bubble emis- 
sion is consistent with emission from a thermal plasma 
in collisional ionization equilibrium with a tempera- 
ture log(TLB/K) = 6.06lg'o4 and an emission measure 
J n^dl = 0.018 cm~ 6 pc. The temperatures of the Local 
Bubble and the Galactic halo components are in good 
agreement with previous ROSAT measurements of other 
direc tions (|Snowden et al.l fl998. 200Q;; lKuntz fc Snowdenl 
I2000h . However, our Local Bubble emission measure 
is 3-10 times larger t han the i?OS^ T-determined value 
(|Snowden et al.lll998l ). This discrepancy is explained by 
the fact that we use a different plasma emiss ion code and 
abundance table from lSnowden et al.l (|1998l ). 

2. Our Local Bubble temperature disagrees with 
the results of XMM-Newton observations of the Local 
Bubb le in other directions, wh ich find log(TT ; R/K) & 
6.2 (jFrevberg fc Breitschwerdtl 120031: iFrevbergj 120041: 



iFrevberg et alll200l . If we use a model with this higher 
Tlb, we find it over-predicts the Galactic halo O VI 
emission by several orders of magnitude. This therefore 
suggests that the Local Bubble is thermally anisotropic. 
However, it is possible that for some of these other XMM- 
Newton observations the foreground emission is being 
contaminated by non-Local Bubble emission from Loop 
I. 

3. Our data are also consistent with a non- 
equilibrium model in which the plasma is underion- 
ized. However, while an overionized recombining plasma 
model is observationally acceptable for certain densi- 
ties and temperatures, it generally gives a very young 
age for the Local Bubble: < 6 x 10 5 yr. This 
is several tim es lower than the Local Bu b ble a ge in 
the models of IBreitschwerdt fc Schmutzlerl (| 19941 ) and 
IBreitschwerdt et al.l (|1996l) . Such a young age is implau- 
sible, as it would require the Local Bubble to have burst 
highly supersonicaUy out of its natal cloud, meaning that 
a purely recombining model is not self-consistent. 

As stated in the Introduction, X-ray spectroscopy is es- 
sential for distinguishing between models of Local Bubble 
formation. The XMM-Newton observations presented 
here have added to the evidence against overionized, re- 
combining models of the Local Bubble. Forthcoming 
Suzaku data from these same viewing directions will ex- 
tend our sensitivity to lower energies, enabling us to place 
constraints on the carbon and nitrogen line emission in 
the ^0.3-0.5 keV energy range. This increased spectral 
information will help us constrain the other major class 
of non-equilibrium Local Bubble model, namely an un- 
derionized plasma that is in the process of ionizing. 
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APPENDIX 

MEASURING T LB FROM FREYBERG ET AL.'S (2004) BARNARD 68 DATA 

IFrevberg et alJ ([2004) do not quote a value for Tlb measured from their XMM-Newton observation of the Bok 
globule Barnard 68. Here we describe how a Local Bubble temperature may be inferred from the data they do present. 

Barnard 68 casts a deep shadow in the soft X-ray background. However, the shadowing is not complete at O vni 
energies, implying at least some of the O VIII em iss ion or iginates in front of Barnard 68. If one attributes this emission 
to the Local Bubble, then from IFrevberg et all s (2004) Figure 2 (which gives the ratio of the on-cloud to off-cloud 
intensity as a function of energy) one can estimate the O vn:0 vni intensity ratio for the Local Bubble emission, and 
hence infer a temperature. 

The ratio Rq viii of on-cloud to off-cloud intensity for O vni emission is 



i?o 



L 



VIII 



O VIII 



1 Bq viii 



-^OVIII + c 



1 Bq vni 



(Al) 



where Loviii and fiovin are the Local Bubble and background O viii intensities , and TQ n VIII and Tovm are the on- 
and off-cloud optical depths at O Viii energies. There is a similar equation for Rq vii, the corresponding ratio for O vn 



9 Note that Bq viii will in fact consist of O VIII emission from the Galactic halo, and also continuum emission from the extragalactic 
background. 
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6.2 6.3 
log (T LB /K) 

Fig. 9. — Hi stogram of 1000 randomly generated values of log(TLB/K), based upon the XMM-Newton observation of Barnard 68 
(Frc ybcrg ct al. 2004). The temperatures were calculated from 1000 values Lq^iu/Lq yu, which in turn were calculated using equa- 
tion l lA2t assuming 20% errors on Bovill/Bovn, t~oViii> an d T ovir 

The on- and off-cloud column densities are 



omission. 



TO 23 and 



dFrevberg et~aI1l2004h . We calculate 



absorp tion cross-section s per hy drogen atom using t he lBaiuciriska-Church fc McCammonl [1992j cross-sections (except 
for He; lYan et aT]|1998[ ) with the I Wilms et~aTl [2000l | interstellar abundances. The values at 0.57 keV (O vn) and 0.654 
keV (O vin) are 6.96 x 10~ 22 and 4.74 x 10~ 22 c m 2 , respectively. Given the large on-cloud optical depths, we can 
ignore the term involving TQ n VIII in equation (|A1[) . and the corresponding term involving Tg n VII in the equation for 
.Rovii- Hence, by combining equation (|A1|) with its -Rqvii analogue we find that 



^OVIII 



-Bqviii-Ro vine" 



off 
r O VIII 



Lo 



VII 



1 — i?o VII 

l--Ro 



VIII 



(A2) 



-Bovii-Roviie 

Figure 2 in lFrevberg et al] (2004) plots the on-off intensity ratio every 50 eV. Thus, for the purposes of this estimate, 
we read off -Rovii and -Rovm at 0.55 and 0.65 keV, respectively (yielding i?ovn = 0.71 and -Rovm = 0.51), and we 
take the O vn and O VIII "bands" to be 0.525-0.575 and 0.625-0.675 keV, respectively. If we assume the Galactic 
halo and extragalactic background are isotropic, we can use our "standard" model results in Table [2] to estimate that 
-Bo vm /-Bo vn = 0.455. We therefore obtain Lo viii/Bovn = 0.20. By using XSPEC to calculate APEC model fluxes 
in the above energy bands, we find that this ratio corresponds to log(TLB/K) = 6.24 (we obtain the same value if 
we use MeKaL model fluxes, instead of APEC). By using a full plasma emission model to calculate flux ratios, as 
opposed to simply using the O vn and O Vin line emissivities, we also take into account the contribution of continuum 
emission to the two energy bands being considered. 

This estimate of Tlb assumes that the value of Bo vin/-Bo vn measured from our XMM-Newton spectra can be 
applied to this pointing direction. We have also assumed a single absorption cross-section for each energy band, when 
in fact the absorption cross-section will vary across each band. To estimate how much of an effect these assumptions 
have on our result, we used a Monte Carlo method similar to that described in $5.2\ to calculate 1000 values of Tlb 
assuming 20% errors on Bqyiii/ Bq vii, T o n viiD an d t ovip Visual inspection of the resulting histogram of values 



indicates that this estimate of Tlb is accurate to within 0.1 dex (see Fig. [5]). 
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